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ABSTRACT 

We analyze the propagation of waves in sunspots from the photosphere to the chromosphere using 
time series of co-spatial Call H intensity spectra (including its line blends) and polarimetric spectra 
of Si I A 10827 and the He I A 10830 multiplet. From the Doppler shifts of these lines we retrieve the 
variation of the velocity along the line-of-sight at several heights. Phase spectra are used to obtain 
the relation between the oscillatory signals. Our analysis reveals standing waves at frequencies lower 
than 4 mHz and a continuous propagation of waves at higher frequencies, which steepen into shocks 
in the chromosphere when approaching the formation height of the Call H core. The observed non- 
linearities are weaker in Call H than in He I lines. Our analysis suggests that the Call H core forms 
at a lower height than the He I A 10830 line: a time delay of about 20 s is measured between the 
Doppler signal detected at both wavelengths. We fit a model of linear slow magnetoacoustic wave 
propagation in a stratified atmosphere with radiative losses according to Newton's cooling law to the 
phase spectra and derive the difference in the formation height of the spectral lines. We show that 
the linear model describes well the wave propagation up to the formation height of Call H, where 
non-linearities start to become very important. 
Subject headings: Sun: magnetic fields; Sun: oscillations 



1. INTRODUCTION 

Since the first dete ction of waves in sunspots 
(Be ckers Tallantl [T969) , many studies have been car- 
ried out to understand the physics of these waves, from 
the observational as well as theoretical point of view. 
Waves in sunspots are different from those observed in 
quiet Sun due to the presence of the magnetic field. 
They show a variety of behavior depending on the height 
and the region of the sunspot where they are observed. 
Usually, these waves are classified as photospheric um- 
bral oscillations, chromosph eric umbral oscillations, and 
running penumbral waves (jLitesI fl992l ). However, all 
these kind of waves seem to be different manifesta- 
tions o f the same global propagation of ma gneto-acoustic 
waves (jRouppe van der Voort etal] |2003). The reasons 
to study waves extend beyond the derivation of their 
properties, because waves can also be used as indepen - 
dent diagnostic of sunspot structure (jDuvall et al.lll996l ) . 

At the photosphere in the umbra, the power spectra 
of the oscillations are quite similar to the corresponding 
ones of the quiet Sun, with a broad distr ibution of fre- 
quen cies and a clear peak at about 3 mHz (jThomas et al.l 
Il982f ). In sunspots, these oscillations in the 5 min 
band have reduced amplit udes compared to th e sur- 
rounding quiet photosphere (jAbdelatif et al.lfl986). The 
magnetic field also modifies the p— modes and produces 
new modes of fluctuations not present in the quiet Sun 
(jCallv fc Bogdarill997t IKhomenko et aTll2009h . Gener- 
ally, most of the p hotospheric umb ra is covered by co- 
herent oscillations (jKobanovlll990[ ). Waves are usually 
studied by measuring the fluctuations of intensity and 
line-of-sight velocity from Doppler shift, but they are also 
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suppo sed to affect magn etic field. From full Stokes inver- 
sions, iLites et all (|1998[ ) found an upper limit of 4 G for 
the amplitude of 5 minute oscillat i ons in magnetic field 
strength, while iBellot Rubio et all (j2000l) detected vari- 
ations around 7-11 G. In a later work, IKhomenko et al.l 
(|2003l ) interpreted magnetic field oscillations as being due 
to fast and slow magneto-hydrodynamic (MHD) wave 
modes, producing both intrinsic variations of the mag- 
netic field and those due to the shift of the line forming 
region. 

The chromosphere of sunspots is dominated by 3 min 
oscillations. The power spectra usually show sharp peaks 
around 5-6 mHz, and the power drops gradually to noise 
value s in the interval from 8 mHz to 10-15 mHz (jLitesI 
119841 ). Several spectral lines, formed at different heights 
from the photosphere to the chromosphere, can be used 
to sample the wave propagation properties. Simultane- 
ous time-series of these spectral lines are a powerful tool 
for studying sunspot waves. From phase spectra between 
the umbral oscillations ob served in the photosphe ric line 
Fel A 5233 and in Ha, iGiovanelli et al.l (|1978l ) found 
that the phase d elay indicates upward wave propaga- 
tion. iLitesI (|1984f ) inferred that slow mode waves propa- 
gate upward into the chromosphere in the frequency band 
around 6.5 mHz, based on the phase differences between 
the oscillations of Call A 8498, Call A 8542 and Fel A 
5434. From the study of the Doppler shifts obser ved in 
the in tensity profiles of the He I A 10830 multiplet. ILitesI 
(|1986l ) pres ented evidence o f shoc k formation at the chro- 
mosphere. iCenteno et al.l (j2006h reproduced the phase 
spectra between chromospheric and photospheric veloc- 
ity oscillations with a model of linear vertical propagation 
of slow magnetoacoustic waves in a stratified magnetized 
atmosphere that accounts for radiative loses, finding a 
good agreement between the theoretically computed time 
delay, and the one obtained from the cross-correlation of 
photospheric Si I A 10827 and chromospheric He I A 10830 
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velocity maps, filtered around the 6 mHz band. They 
showed that the chromospheric 6 mHz signal is a result 
of linear wave propagation of the photospheric perturba- 
tions in the 6 mHz range, rather than the consequence of 
the nonli near interaction of ph o tosph eric modes as pro- 
posed bv iGurman fc Leibacherl (j!984f ). 

The works cited above were limited to the study of os- 
cillations at only two heights (one photospheric and one 
chromospheric), separated by several hundreds of kilo- 
meters. It is thus interesting to perform a more detailed 
sampling of the sunspot atmosphere, using more spec- 
tral lines formed at several intermediate heights between 
these two regions. On the one hand, observationally 
detected spatial wave patterns in sunspots are rather 
complex (Bogdan & Judge 2006). On the other hand, 
recent numerical simulations of waves in sunspots also 
suggest a complex picture of the fast and slow magneto- 
acoustic waves propagating simultaneously in the same 
layers b ut in different directions and w ith different phase 
speeds (|Khomenko fc Colladosl 120091 ) . This requires a 
more refined multi-layer study of sunspot waves. Studies 
of this kind often represent an observational challenge 
since several spectral lines have to be detected simulta- 
neously not only in intensity but also in polarized light. 
In our paper, we report on such multi-line spectropolari- 
metric observations. Our aim is to cover the gap between 
the photospheric and chromospheric signals and analyze 
sunspot oscillations at the formation heights of several 
spectral lines formed between Si I and He I. For that we 
use simultaneous observations obtained with two instru- 
ments, the POlar imetric LIttrow Spectrograph (POLIS, 
iBeck et a l. 2005b ) and the Tenerife Infrared Polarimeter 
II (TIP-II, iCollados et all 120071 ). attached to the Ger- 
man Vacuum Tower telescope at the Observatorio del 
Teide at Tenerife. Apart from information about waves, 
our multi-layer study has also allowed us to estimate the 
formation heights of the spectral lines in the sunspot at- 
mosphere, including Call H line, several Fel blends in 
the wing of this line and the infrared lines of Si I at A 
10827 A and He I at A 10830 A. 

The structure of the paper is the following. In Sect. [2j 
the observations and data reduction are explained. Sec- 
tion [3] describes the analysis of the velocity oscillations 
at several heights. The results are discussed in Sect. HJ 
which also presents our conclusions. 

2. OBSERVATION AND DATA REDUCTION 

The observations analyzed in this work were obtained 
on 2007 August 28 with two different instruments, PO- 
LIS and TIP-II, attached to the German Vacuum Tower 
Telescope (VTT) at the Observatorio del Teide. Simul- 
taneous and co-spatial scans of a sunspot located near 
the center of the Sun (/i = 0.96) were taken with both 
instruments. The slit was placed over the center of the 
sunspot. The observations were obtained with real-time 
seeing correcti on by the Kiepenheuer-Ins titute adaptive 
optics system (jvon der Luehe et al.ll2003l ). 

The spectra of the blue channel of POLIS include 
Call HA 3968 A intensity profiles and some photospheric 
line blends in the wings of Call H, covering a spectral 
range from 3964.9 A to 3971.3 A with a spectral sam- 
pling of 20 mA pixel -1 and a spatial sampling of 0"29 
per pixel. The Ca spectra were reduced for the flat- 
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Fig. 1. — Intensity of IR continuum. Vertical lines mark the 
position of quiet Sun-penumbra and penumbra-umbra boundaries. 
Time increases from bottom to top. 

field (B eck et al.l l2005ab). and were also corrected for 
the transmission curve of the order-selecting interference 
filter in front of the camera. For the wavelength calibra- 
tion, the line-core positions of the iron lines at 3965.45, 
3966.07, 3966.63, 3967.42 and 3969.26 A in an average 
quiet Sun region were determined by a second order poly- 
nomial fit. We then determined the wavelength scale that 
matched best all the position values. 

TIP-II yielded the four Stokes parameters IQUV in a 
spectral range from 10822.7 A to 10833.7 A with a spec- 
tral sampling of 11 mA pixel -1 and a spatial sampling 
of 0"18 per pixel. This spectral region contains informa- 
tion about two different heights of the solar atmosphere 
due to the presence of two spectral lines. The Si I line 
at 10827.09 A is formed at photospheric heights, whereas 
the He I A 10830 triplet, which includes a weak blue com- 
ponent at 10829.09 A and two blended red components 
at 10830.25 and 1083 0.34 A, forms in the chromosphere 
(jCenteno et al.l 1200 8) . In this case, the wavelength cal- 
ibration was done using the Si I and He I lines as refer- 
ences. 

The full data set consists of three temporal series; each 
of them lasting about one hour. Three scan steps with 
0"5 step width were taken for the two first series, while in 
the last series only two spatial positions were used. The 
cadence was different for all series. Table [1] shows the 
time step between two spectra taken at the same spatial 
position, the time when the observation were obtained 
and the number of repetitions of each series. 



TABLE 1 

Summary of observations: 





tstart (UT) 


tend (UT) 


At (s) 


N spectra 


Series 1 


07:55:02 


08:54:48 


21 


170 


Series 2 


09:01:02 


10:15:25 


18 


250 


Series 3 


10:43:34 


11:58:43 


7.5 


600 



Due to the differential refraction in the earth atmo- 
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Fig. 2. — Temporal evolution at one position in the umbral region of series #2. Left: Si I and He I intensity; center: Si I and He I Stokes 
V; right: Call H intensity. The horizontal axis represents wavelength, with the origin at the position of the Si I A 10827 A rest wavelength 
in (a) and (b) and at the position of the Call H A 3968 A rest wavelength in (c). The vertical axis represents time, increasing from bottom 
to top. 

in the infrared continuum. The regions of quiet Sun, 
penumbra, and umbra of the sunspot are well defined, 
the vertical lines indicate the boundaries between these 
areas. The temporal evolution of the Stokes I and Stokes 

V spectra for the TIP data and the intensity around the 
Call H core in the POLIS data are plotted in Fig. [2] at 
a fixed position inside the umbra. The Stokes / and V 
profiles from TIP contain the Si I line at AA =0 A (the 
rest wavelength of the silicon line was determined from 
the quiet Sun region and was set as the origin) and the 
He I line at AA =3.25 A. The helium line profile shows 
periodic shifts with large displacements to the blue and 
the red. The core of the Call H line shows a strong 
emission peak inside the umbra. The Doppler shift of this 
emission peak dev elop a saw-tooth pattern, similar to the 
He I line (see also lRouppe van der Voort et al .1120031 ). 

3. DATA ANALYSIS AND RESULTS 

In this paper, we focus on the line-of-sight (LOS) veloc- 
ities. For all the spectral lines besides Call H, Doppler 
velocities were inferred by measuring the position of the 
intensity minimum. The wavelengths close to the core of 
the line were fitted with a second order polynomial. The 
location of the minimum of the parabola was taken as 
the line-core position. This procedure was performed to 
obtain the Doppler shifts of Si I, He I, and the Fel lines 
at 3965.45, 3966.07, 3966.63, 3967.42 and 3969.26 A. In 
the case of the TIP data, the Doppler shift of the Stokes 

V zero crossing, where the polarization signal intersects 
the zero level, was derived as well. The Doppler shifts 
from the intensity and Stokes V profiles are very similar 
due to the large magnetic filling factor in the sunspot 
umbra. 

The behavior of the Call H core is different. It exhibits 
a prominent peak at the center of the line in highly mag- 
netized regions (top line of Fig. |4]), while in a field- free 
atmosphere the center of the line has a minimum between 
two lobes at both sides with their co rresponding maxim a 
(bottom line of Fig. |4]) (see also, e.g. JLiu & Smithlll972f ). 
In the umbra, the Doppler shift was retrieved from the 
spectral displacement of the maximum of the core emis- 
sion. In the low magnetized, region it was obtained from 
the shift of the central minimum. All in all, we have 
obtained maps of the LOS velocity at each spatial point 



Fig. 3. — Temporal variation of the theoretical spatial displace- 
ment perpendicular to the slit due to differential refraction (Solid 
line). Asterisks mark the values measured at the telescope. 



sphere (e.g. JFilipp enko 1982), the spectra of POLIS and 
TIP-II are not fully co-spatial. The spatial displacement 
of the two wavelengths (3968 A, 10830 A) perpendicular 
to the slit depends on the date and time of the obser- 
vations, the slit orientation, and the l ocation of the firs t 
coelostat mirror (see Appendix A of Beck et al.l 12008). 
On the first day of the observation campaign, we took 
a set of large-area scans at different times for an accu- 
rate determination of the displacement. The solid line in 
Fig. [3] shows the theoretically predicted spatial displace- 
ment perpendicular to the slit due to differential refrac- 
tion, while the asterisks are the measured displacements; 
the match between both is remarkable. To guarantee an 
overlap between the observations in the two wavelengths, 
we thus positioned the scan mirror inside of POLIS at the 
beginning of each observation such that it compensated 
the spatial displacement for about the middle of the ob- 
servation. Moreover, small repeated scans of 2-3 slit spa- 
tial positions, separated by 0"5, were taken in order to 
sample a wider region of the Sun and prevent possible 
errors between the theoretical differential refraction and 
the actual one. 
Figure [T] shows an intensity map from a wavelength 
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Fig. 4. — Averaged intensity profile of the Call H line at different 
locations. From bottom to top: quiet Sun, penumbra, and umbra. 
The intensity scale is different for all of them. The profiles have 
been displaced from each other. The Fel lines analyzed in this 
work are marked. 

covered by the slit at the formation heights of the eight 
spectral lines, except for He I, whose line depth in the 
non-magnetized regions is too small to determine its lo- 
cation. 

We have selected the scan steps which give the best 
alignment of the data from both instruments to obtain 
a better correction of the differential refraction. The use 
of different steps introduce a systematic delay between 
the oscillatory signals from POLIS and TIP-II, since the 
same spatial location is observed with a time lag. How- 
ever, in the case of the data from series 3, the theoretical 
correction of the differential refraction worked perfectly 
and the optimal alignment was found using the simul- 
taneous scan steps. This data set has been chose n to 
calculate the phase difference spectra from Section 13.31 
An additional finer alignment along the slit between TIP 
and POLIS was done by a cross-correlation of the veloc- 
ity maps of two spectral lines formed at similar heights, 
using the line pair Call H and He I for the chromosphere, 
and Si I and Fe I A 3969.26 A for the photosphere. Finally, 
we resampled all velocity maps with the sampling of the 
POLIS Call H data, i.e., 0"29 per pixel. 

3.1. LOS velocity 

Figure [5] shows the temporal evolution of the LOS ve- 
locities in the sunspot region obtained from the Doppler 
shift of several spectral lines, sorted by formation height. 
Negative velocities (appearing as black shaded regions) 
indicate upflows, where matter is approaching the ob- 
server, while white regions are downflows. Figures [5ja-d) 
reveal a similar pattern since all these lines are formed at 
photospheric heights. The velocity maps of the Fel line 
blends of Call H have lower quality than the Si I maps, 
especially inside the umbra, due to the coarser spectral 
sampling and the small line depth of the lines. 

The velocity field at chromospheric heights is given in 
Fig. [St (Call H) and[5]F (Hel). Both figures show a sim- 
ilar wave pattern and similar amplitudes, but, as a more 
detailed analysis will reveal, the velocity oscillations of 
the He I line are slightly retarded relative to those in Ca. 
The wave pattern measured at the chromosphere in the 
umbral region has a smaller spatial coherence than that 
observed in the photosphere. It also differs in its larger 



peak-to-peak variations of about 8 km s _1 and its period 
of about three minutes. 

In the region of the penumbra, which can be seen 
in Fig. [Sf at a position between 19" and 23" and be- 
tween 34" and 39", there is also a characteristic pat- 
tern of alternating positive and negative velocities, but 
with lower amplitudes than the waves in the umbra and 
with longer periods, corresponding to running penum- 
bral waves (see, e.g.JGiovane 

iBloomfield et al.ll2QQ7| ). Thev start at the inner penum- 
bra, and their wave front is delayed in the regions of 
the outer penumbra, so their propagation through the 
penumbra appears as a slope in the diagram. An exam- 
ple of such a slope can be seen in an area between 19" 
and 23" during the first 22 minutes (white rectangle). 
This slope is smaller near the outer penumbra, indicating 
a decrease of the propagation speed as the disturbance 
travels from the inner to the outer penumbral bound- 
ary. The estimation of the velocity for the Call H line 
core (Fig. [5^) in the penumbral region is poorer due to 
the change of the shape of the line that does not have a 
prominent emission peak any more (see middle curve in 

Fig. mi. 

The left panels of Fig. [6] show the histograms of the 
LOS velocities in the umbra of the sunspot obtained from 
all the spectral lines used in this study. We fitted a 
Gaussian to the distributions to estimate the root-mean- 
square (rms) velocity value. The histograms are sorted 
from bottom to top with increasing rms velocity, with Si I 
having the smallest and He I the largest rms velocity. The 
Fel lines at 3966.0 and 3966.6 A have an identical rms 
velocity value, as they are formed at close heights. This 
is also suggested by their nearly identical line depth (see 
Fig. Hj). From these histograms and the variation of the 
rms velocities, we obtain the first estimate of the relative 
formation height of the spectral lines since we expect a 
monotonic increase of the rms velocity with height. 

For comparison, the right panels of Fig. [6]show the his- 
tograms of LOS velocities for the same lines in a region of 
the quiet Sun. The He I histogram is not shown because 
the He I absorption in quiet Sun is too low to retrieve a 
velocity. The quiet Sun histograms also show the increase 
of the rms velocity with the formation height, although 
the order of the rms velocity slightly differs from the 
one obtained in the sunspot, since all the Fel lines form 
in a thin layer. As the oscill atory power in sunsp ots is 
suppressed (see, for example. iHindman et al.l ll997). the 
rms velocities in the quiet Sun are larger. Another pos- 
sible source of the broadening of the quiet Sun velocity 
histogram can come from the granulation velocities not 
present inside the umbra. We did not remove the gran- 
ular component from the velocity variations in the quiet 
Sun and it can contribute to some extent to the overall 
rms v elocity variations (see, e.g. iKostvk fe Khomenkol 

dool). 

Figure shows an example of the temporal evolution of 
the LOS velocity obtained from the chromospheric Hel 
line and the Call H line core, the four photospheric Fel 
lines, and the photospheric Si I line at one location in- 
side the umbra of the sunspot, indicated in all panels of 
Fig. [5] with a vertical dashed line at x = 28". The plots 
are sorted from bottom to top with increasing formation 
height, as retrieved from Fig. [6] A comparison between 
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Fig. 5. — Velocity maps in the sunspot during series #2 for the different spectral lines. The horizontal axis represents the position along 
the slit and the vertical axis represents time. Black color indicates negative velocities (upflow), white color positive velocities (downflow). 
Top row, left to right: Si I A 10827 (a), Fei A 3966.1 (b), Fei A 3966.6 (c). Bottom row, left to right: Fei A 3969.3 (d), Can H 3968.5 (e), 
He I A 10830 (f). Vertical solid lines represents the limits of the umbra. The vertical dashed line denotes the location of the spectra shown 
in Fig. [2 The white rectangle in (f) marks an area of running penumbral waves. 



the bottom and top panels (which show the Si I and He I 
velocities, respectively) reveals the differences in period 
and amplitude of the waves at photospheric and chro- 
mospheric heights. Between these two layers, the rest of 
the spectral lines sample different heights of the atmo- 
sphere. At photospheric heights (panels c-g), the pattern 
of waves is similar and the signals of the lines with higher 
formation heights are slightly delayed (see e.g. the peak 
at t = 49 min marked with a dashed line). We can see 
that higher layers have larger amplitudes (note that the 
top two panels have a different scale for the velocity). 
The temporal evolution of the LOS velocity of the two 
chromospheric lines (panels a-b) is almost identical, but 
the amplitudes are higher in the case of the He I line. For 
instance from minute 10 to 25, it is clearly seen that the 
oscillations measured with the He I line and Ca II H core 
have a saw-tooth profile that indicates the presence of 
a shock wave train, with a slow increase of the velocity 
followed by a sudden decrease. There is a phase differ- 
ence of about 20 s between the Ca II H line core and the 
He I line velocities, in the sense that the oscillatory signal 
reaches the formation height of the Ca II H line core 20 s 
before that of the He I line. 

Table [2] lists the rms velocities in the quiet Sun (ex- 
cept Hel) and in the sunspot for all the spectral lines, 
together with the ratios between the quiet Sun and the 
sunspot velocities. The relative increase of the rms ve- 
locity of the quiet Sun to the rms velocity of the sunspot 
decreases with the height in the atmosphere, from 2.80 
at the photospheric height of the formation of Si I to 1.29 
at the chromospheric height of the formation of Call H 
line core. The rate of the increase of the rms velocity in 
the sunspot is faster, and the difference between veloc- 
ities from lines formed at higher layers is smaller. The 



growth of the amplitude of the oscillations with height is 
scaled by the pressure scale height, Ho = RT/[ig, where 
T is the temperature, R is the gas constant, g is the 
gravity and ji is the mean weight of the atoms. In the 
umbra the temperature is lower and Hq is smaller than 
in the quiet Sun atmosphere, so the amplitude of umbral 
oscillations rises faster. 

3.2. Power Spectra 

Figure [8] shows the normalized average power spectra of 
LOS velocity of the two chromospheric (He I and Call H) 
and two of the photospheric (Si I and Fei A 3969.3) lines 
inside the umbra of the sunspot. We chose this iron line 
because it has better signal to noise and its formation 
height is distant from the layer where Si I is formed. In 
the photosphere (bottom panel), the power is concen- 
trated between 2 and 4 mHz, corresponding to the 5 
minute band, with a maximum peak at 3.5 mHz. Both 
spectral lines peak at the same frequency, although the 
power of the Fei line is slightly higher. The increase of 
the power at frequencies above 4.5 mHz is more impor- 



TABLE 2 

Rms velocities in ms _1 in the sunspot and the quiet Sun, 
and their ratio (3rd column). 



line 


umbra 


QS 


ratio QS/umbra 


Si i 


90 


250 


2.80 


Fei 3965.4 


170 


256 


1.51 


Fei 3966.1 


174 


275 


1.58 


Fei 3966.6 


180 


273 


1.52 


Fei 3969.3 


201 


306 


1.53 


Can H 


510 


659 


1.29 


He I 


726 
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Fig. 6. — Histograms of LOS velocity measured with several spec- 
tral lines. The left column corresponds to the umbra, the right 
column to the quiet Sun. From bottom to top: Si I, Fel A 3965.4, 
Fei A 3966.1, Fel A 3966.6, Fel A 3969.3, Can H and Hei. Solid 
lines represent the Gaussian fit; its width is indicated in each plot. 



tant than the one for frequencies below this value. 

The velocity power spectra of both chromospheric lines 
(top panel) have a broad distribution of frequencies, with 
the largest power being in the band from 5 to 10 mHz. 
The chromospheric power spectrum has a maximum at 
6.2 mHz an d several se condary peaks around it (see for 
comparison iLitesI Il986[ ) . These frequency peaks corre- 
spond to the chromospheric 3 minutes oscillations. At 
the highest peak of the power spectra, both Call H and 
He I have almost the same power, but for those frequen- 
cies with lower power, the power of the He I is increased 
comparing to the Call H. Note that at heights sampled 
by our spectral lines we do not find a continuous tran- 
sition from the peak at 3.5 mHz to the one at 6.2 mHz, 
but rather a discontinuous behavior between the photo- 
spheric and chromospheric power spectra. However, the 
prominent secondary peak around 5.5 mHz in the power 
spectra of the Fel line is much more obvious than the 
corresponding in the Si I power spectra, which could in- 
dicate some transition towards higher frequencies in the 
power spectra as the waves propagate upward from the 
formation height of the Si I line to the Fe I lines. 




3.3. Phase spectra 



Fig. 7. — Velocity at a fixed position in the umbra for different 
spectral lines, sorted by formation height: He I A 10830 (a), Call H 
3968.5 (b), Fei A 3969.3 (c), Fei A 3966.6 (d), Fei A 3966.1 (e),Fei 
A 3965.4 (f), Si I A 10827.1 (g). The vertical dashed line at t = 49 
min marks a prominent velocity peak. 

A phase diagram gives the phase difference (A0) be- 
tween two signals. In our study, we use to measure 
the time delay between the oscillatory velocity signals 
from two spectral lines and assume that the difference 
between them is mainly due to the difference of the for- 
mation height of the two lines. In the following, we will 
show the phase difference spectra between different com- 
binations of pairs of spectral lines used in this work. To 
obtain the phase spectra, we treated each spatial point 
separately and calculated the Fourier-transform of the 
temporal evolution of the respective velocities. We de- 
rived the phases, and from them the phase difference of 
the two signals as a function of the frequency. There is 
a 2tt ambiguity in the computation of the phase value, 
so all phase differences have been projected in the range 
d=7r. Then we calculated histograms of the relative oc- 
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Fig. 8. — Average umbral power spectra of the LOS velocities. 
Top: chromospheric lines, He I (solid line) and Call H (dashed 
line), both are normalized to the maximum power of the He I line; 
bottom: photospheric lines, Si I line (solid line) and Fel 3969.3 A 
(dashed line), both are normalized to the maximum power of the 
Fel line. 

currence of a given value of the phase differences at each 
frequency taking int o account all the co rresponding spa- 
tial points (see also Krii ger et all I2QQ1L and references 
therein). We obtained the data displayed in Figs. I9HT51 
In addition to the phase difference spectra, we calcu- 
lated the coherence spectra. They provide an estimate 
of the statistical validity of the phase and power spec- 
tra. Considering n pairs of signals Xk(t) and yk(t), whose 
Fourier transforms are Xk(w) and Yk(<u), respectively, 
the coherence is defined as 



Pxy{u) 



(i) 



where A(/)k(uj) = <j> x k(w) — <t*yk(w). In our case, the sub- 
index k covers the spatial position. The coherence eval- 
uates statistically for every frequency lj the relation of 
the A0/ c (cj) for the n ^-signals. It takes the value 1 when 
Acj)k(uj) is the same for all the k. If the phase difference 
of the different k is arbitrary, the coherence takes very 
low values. We selected a confidence limit at 0.7, and for 
frequencies with coherence above this value we consider 
the phase spectra to be reliable. 

We also analyzed the increase of the amplitude of the 
oscillations. We calculated the amplification spectra as 
the ratio between the power at two layers, both of them 
averaged all over the umbra: 



A 



(2) 



3.3.1. Theoretical model 

Following Ce nteno et al.l (12006V the observations were 
compared with a model of linear vertical propagation of 
slow magneto-acoustic wave in an isothermal atmosphere 
that includes radiative losses described by Newton's cool- 
ing law. Assuming that the amplitude of the vertical 
velocity changes with height by 



V(z) = V e 
the dispersion relation for such waves is 



2 ~ 2 



(3) 



(4) 



where 



u ac = c /2Hq,c = <ygH ,<y = - : , (5) 

1 - IUTr 

and tr is the radiative cool ing; time for an optically thin 
perturbation (jSpi egel 1953): 



tr = P c v /(16 X (TrT 3 ). 



(6) 



The phase difference between oscillations at two 
heights is calculated as the difference of the imaginary 
part of the argument of the complex exponential in Eq.[3J 
that is, Acf) = kRAz, where Az is the geometric dis- 
tance between the two heights and kR is the real part of 
k z . The amplification spectrum is given by the ratio be- 
tween the amplitude at the two layers and is obtained as 
A = e 2 ( 1 /( 2 ^o)-fe/)A^^ with k\ being the imaginary part 
of k z . 

This model allows us to fit the phase and amplification 
spectra with three free parameters: the temperature of 
the atmosphere, T, the difference in height between two 
lines, Az, and the typical time scale in which the tem- 
perature fluctuations are damped radiatively, tr. These 
free parameters are manually tuned to match the effec- 
tive cut-off frequency and the slope (including its vari- 
ations) of the phase difference spectra in the regime of 
propagating waves above the cut-off frequency. 

3.3.2. Sil-Hel phase spectra 

Figure [9] shows the phase difference between the ve- 
locity signals measured in the photospheric Si I line and 
the chromospheric He I line. The phase difference is zero 
for frequencies below 4 mHz. At these frequencies the 
coherence is high, except in the range between 1 and 2.5 
mHz, and the atmosphere oscillates as a whole, i.e., the 
waves are stationary. From 4 mHz to 7 mHz, the phase 
difference increases linearly with the frequency and the 
coherence is (more or less) above the confidence limit. It 
indicates that waves at these frequencies propagate from 
the photospheric layer, where the Si I line forms, to the 
chromospheric layer, where He I forms. The phase differ- 
ence spectra for higher frequencies are very noisy, and no 
meaningful conclusions are possible. 

The bottom panel of Fig. [9] shows the ratio of chromo- 
spheric to photospheric power as a function of frequency. 
We fit both the phase difference and the amplification 
spectra simultaneously. The parameters of the fit are 
listed in Table OH The solid red line in the phase diagram 
in Fig. ([9]) represents the ph ase dif ference according to 
the model presented in Sect. 13.3. 1\ calculated with the 
parameters that best fit the observations. The theoreti- 
cal amplification spectra matches rather well to the ob- 
servational one, with an order of magnitude agreement 
in the amplification factor. 
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Fig. 9. — Top: Phase spectra between the LOS velocities of the 
Si I and the He I lines. The color code shows the relative occur- 
rence of a given phase shift. The red line represents the best fit 
from the theoretical model. Middle: Coherence spectra. The hor- 
izontal dashed line at 0.7 marks the confidence limit. Bottom: 
Amplification spectra. The red line represents the best fit from 
the theoretical model. 



TABLE 3 

Best-fit parameters of the theoretical model 



Line pair 


T (K) 


Az (km) 


tr (s) 


Si-He 


4500 


900 


45 


Si-Ca 


4000 


650 


45 


Fe 3969.3-He 


4000 


450 


45 


Fe 3969.3-Si 


4000 


-280 


30 


Fe 3969. 3-Si (quiet sun) 


4500 


-200 


30 


Fe 3965.4-Fe 3969.3 


3500 


30 


45 


Fe 3966.0-Fe 3969.3 


3500 


20 


45 


Fe 3965.4-Fe 3966.0 


3500 


10 


45 


Ca-He 


6000 


100 


15 



3.3.3. Sil- C all H phase spectra 

Figure [10] shows the phase, coherence and amplifica- 
tion spectra between the velocity measured with the Si I 
line and the core of Call H line. Frequencies below 1 
mHz have a very low coherence, so these values are not 
reliable. Surprisingly, frequencies in the range between 
1 and 2.5 mHz present an anomalous behavior, with a 
phase difference around 2.5 rad and a high amplifica- 
tion, while their coherence is remarkably high. The long 
period of waves with these frequencies hinders their anal- 
ysis, so further observations with longer temporal series 
are required to study this behavior. The rest of the phase 
spectra is similar to the one between Si I and He I, with 
zero phase difference between 2.5 and 4 mHz, indicating 
stationary waves, and an almost linear increase of the 
phase difference between 4 and 7 mHz, corresponding to 
upwards propagating waves. In all this frequency range, 
the coherence is high. However, in this case the slope of 
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Fig. 10. — Phase, coherence and amplification spectra between 
Si I and Call H. The format of the figure is the same as Fig. [9] 

the increase of the phase difference is smaller than in the 
previous one, and the amplification also has lower values. 
It means that the phase delay between the velocities at 
the formation height of these two lines is smaller than 
the delay between Si I and He I and the amplitude of the 
waves at the formation height of He I is larger than that 
of waves at the layer where Ca II H forms. Both the phase 
and amplification spectra locate the formation height of 
the Call H core below the He I line. 

The parameters retrieved from the fit of the phase and 
amplification spectra to the theoretical model are listed 
in Table 02 Comparing the ones retrieved for the line 
pair Sil-Hel with the pair Sil-Call H, the Call H core 
formation height in the umbra is around 250 km below 
the He I line. The temperature obtained for the phase 
difference between Sil and Call H is smaller by 500 K; 
waves traveling from the formation height of Ca II H core 
to that of He I pass presumably through a region where 
the temperature is increasing. 

3.3.4. Fei phase spectra 

To study the properties of oscillations at photospheric 
heights, phase diagrams between pairs of the photo- 
spheric lines were calculated. We can assume that the 
Fei lines in the Call H wing form at three different 
heights according to their line depth (Fig. 3) and to the 
width of the velocity histograms (Fig. [6]). We take the 
lines Fei A 3965.4, Fei A 3966.0 and Fei A 3969.3 as rep- 
resentative of these heights, since Fe I A 3966.6 and Fe I A 
3967.4 seem to form at a similar height as Fei A 3966.0. 
As an example, Fig. fTTI shows the phase spectra, coher- 
ence and amplification spectra obtained between between 
Fei A 3966.0 and Fei A 3969.3. The phase spectrum 
shows that the phase difference is almost zero for all the 
frequencies that can be trusted according to the coher- 
ence spectra (from to 10 mHz). We conclude that we 
can not retrieve the phase shift between the Fe I lines re- 
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Fig. 11. — Phase, coherence and amplification spectra between 
Fei A 3966.0 and Fei A 3969.3. The format of the figure is the 
same as Fig. [9] 

liably, as the geometrical difference in their heights of the 
formation is too small. However, the amplification spec- 
tra reflects some increase of the amplitude with height. 
We managed to fit both the phase and the amplifica- 
tion spectra with the wave model described in the previ- 
ous section. These fits yield, indeed, a small geometrical 
height difference between the Fe I lines. The temperature 
and the cooling time obtained from the fit are identical 
for all pairs of iron lines (Table [3]). 

3.3.5. Fel-Sil phase spectra 

The phase difference spectrum between the Fei A 
3969.3 line and the Si I A 10827 line is shown in the top 
panel of Fig. [12j For frequencies below 2 mHz, the phase 
spectrum is very noisy and has no coherence (see middle 
panel) , indicating that waves at the heights of formation 
of these two lines are not related. At frequencies in the 
band of 2.5 — 4 mHz, some oscillatory power is present, 
the coherence is high and the phase difference is close to 
zero, so these waves are evanescent. The phase difference 
spectrum for high frequencies in the range from 4 to 9 
mHz shows a decreasing tendency, indicating that the 
waves reach the formation height of the Si I line before 
the Fei line. This phase spectrum was also fitted with 
the wave propagation model. The solid red line in top 
and bottom panels of Fig. fT2l shows the result of this fit. 
The amplification spectrum is not reliable at frequencies 
below 3 mHz due to the low coherence, but above this 
value the agreement of the theoretical and observational 
spectra is good, indicating that the amplitude of the os- 
cillations measured in Si I is about twice lower than the 
amplitude measured in the Fei 3969.3 line. 

As a summary of the propagation at photospheric 
heights, we conclude that all observed photospheric spec- 
tral lines (five Fei lines and the Si I line) fluctuate with 
a dominant period of 5 minutes. While all the Fe I lines 
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Fig. 12. — Phase, coherence and amplification spectra between 
Fei 3969.3 and Sii. The format of the figure is the same as Fig. [9] 

have similar velocity amplitudes, with a maximum peak- 
to-peak amplitude of around 800 m s _1 , the velocity 
amplitude measured in Si I is clearly smaller (maximum 
peak-to-peak amplitude around 400 m s _1 ). According 
to the phase diagram between the Fei lines (Fig. fTTj) 
and the amplification and phase diagram between a Fe I 
line and Si I (Fig. [12]), we conclude that waves with 
frequencies above 4 mHz propagate upwards at photo- 
spheric heights. They first reach the height where Si I is 
formed. As they propagate upward, their amplitudes in- 
crease due to the density fall-off. Then the waves reach 
the formation height of all the Fei spectral lines, that 
all come from a thin layer and for this reason all show 
similar amplitudes. 

3.3.6. Fei-Hel phase spectra 

According to the results extracted from the velocity 
statistics and the phase spectra presented before, the 
Fei lines form at some height between the formation 
height of the Si I line and the chromospheric He I line and 
the Call H core. This means that they give information 
about a high photospheric layer, located at an intermedi- 
ate height in the propagation of the waves from the pho- 
tosphere to the chromosphere. Figure [T3l shows the anal- 
ysis of the phase differences and amplification between 
the velocity signal measured at the formation height of 
the Fei A 3969.3 line and the He I line. The phase spec- 
trum is similar to the other two between photospheric 
and chromospheric lines (Si I- He I and Sil-Call H). For 
frequencies between 2 and 4 mHz, the phase difference is 
about zero, and it increases with the frequency between 
4 and 7 mHz, but the slope is smaller than for Sil-Hel. 
In the range between 2 and 7 mHz, the coherence is sig- 
nificant. The amplification spectrum also shows a lower 
amplification, compared to the Sil-Hel case. The tem- 
perature, height difference and radiative damping time 
retrieved from the fit are listed in Table [3l 
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Fig. 13. — Phase, coherence and amplification spectra between 
Fel A 3969.3 and Hel. The format of the figure is the same as Fig. 

m 

3.3.7. Call H -Hel phase spectra 

At chromospheric heights, we have the Doppler veloc- 
ities obtained from the cores of the Call H line and the 
Hel A 10830 line. The top panel of Fig. [Ml shows the 
phase difference diagram between them. At frequencies 
between 1-2 mHz, the coherence (middle panel of Fig. 
[H]) is low and the phase spectrum is noisy. The coher- 
ence spectrum shows that the phase spectrum is reliable 
between 2 and 12 mHz. In the frequency band between 
2 and 4 mHz, the phase difference is about 0, indicating 
that there is no propagation and the waves are evanes- 
cent. From z/=4 mHz to z/=ll mHz, the phase difference 
increases, starting from A<fi = and showing a small pos- 
itive slope. It means that as waves propagate upwards, 
they reach the Call H core formation height just before 
the Hel one. The same conclusion was obtained pre- 
viously from comparing the temporal variations of the 
Call H core and Hel velocities presented in Fig. Efa-b), 
where the oscillatory signal of He I is delayed by 20 s with 
respect to the Call H one. 

The differences in the amplitudes between Ca II H core 
and Hel velocities are in line with their phase spectra 
(top panel in Fig. fl4|) . since the amplitude of the Hel 
velocity is bigger than the Call H. A more detailed in- 
spection of the amplification spectra between both chro- 
mospheric lines (bottom panel of Fig. fl4|) reveals that 
the oscillatory signal is amplified between 2 and 4 mHz, 
but there is no amplification for frequencies in the range 
from 4.5 to 10 mHz (the ratio between the amplitudes is 
around unity). Waves in the 2-4 mHz frequency range 
are evanescent, and the increase of their amplitude from 
the photosphere to the chromosphere is not so high, keep- 
ing them in a linear regime. At the high layers, their 
amplitude still increases due to the drop of the density 
with height. On the other hand, waves with frequencies 
between 4 and 10 mHz propagate upwards to the chromo- 
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Fig. 14. — Phase, coherence and amplification spectra between 
Call H and Hel. The format of the figure is the same as Fig. [9] 

sphere and develop into shocks (top panels of Fig. [7]). In 
this non-linear regime, their amplitudes do not increase 
with height, explaining the observational amplification 
spectrum around unity for frequencies in the range 4- 
10 mHz. The model of linear wave propagation in an 
isothermal stratified atmosphere predicts a higher am- 
plification. In the case of waves with frequencies above 
10 mHz, they do propagate upwards, but their ampli- 
tude at the photosphere is so low that they do not reach 
a non-linear regime and their amplitude still increases 
with height at chromospheric layers. All in all, the model 
is not suitable for the description of waves at heights of 
Call H core and He I formation where the non-linearities 
start to become important. 

It is expected that the propagation of non-linear waves 
happens at a higher speed compared to the linear case. 
This would decrease the phase difference between two 
layers. Therefore, we can expect that the height differ- 
ence between the formation layers of Call H and Hel 
lines presented in Table [3] is underestimated by our lin- 
ear model; the value of about 100 km is the lower limit 
of this difference. One of the issues that arises from this 
fact is the evaluation of the height range where the linear 
model of wave propagation is valid. It is clear that it fails 
between the formation height of the Ca II H core and the 
He I line due to the non-linearities. However, in previous 
sections we have applied successfully this model to the 
propagation between the photospheric Si I line and the 
chromospheric Call H core and Hel line, and between 
one of the Fel lines and the Hel line. It means that 
in most of their way from the photosphere to the chro- 
mosphere, the linear regime is a good approximation for 
these waves. At some layer between the formation height 
of the Fel 3969.3 and the Hel lines t he wave propaga- 
tion d eparts from the linear regime (cf. iCarlsson fc Stein] 
H997h . 




3.3.8. Fei-Sil phase spectra in quiet Sun 

We calculated the phase difference spectra for a re- 
gion of quiet Sun around the sunspot as well. In this 
region, the Stokes parameters QUV are below the level 
of noise, so we suppose there is no significant magnetic 
field. Figure [15] (top) shows the phase spectrum between 
the photospheric Fel A 3969.3 and Si I A 10827 lines in 
the quiet Sun. The coherence spectrum shows high val- 
ues for frequencies below 9 mHz. For frequencies below 
2 mHz, the phase diagram is very noisy due to the low 
oscillatory power of the two velocity signals, and we do 
not find any clear relation between them. In the spectral 
range from 2 mHz to 4 mHz, the phase difference takes a 
constant value around zero, while for frequencies above 4 
mHz, Acf) decreases with increasing frequency. However, 
the low power at the photosphere at frequencies between 
6 and 8 mHz (Fig. [8]) results in a noisy phase spectrum 
in this band. 

Comparing the phase difference diagrams between the 
same two lines inside the umbra of the sunspot (top panel 
in Fig. [T2|) and in the quiet Sun, we find that in quiet Sun 
the negative slope of the phase difference is less steep. We 
also fitted the quiet Sun phase difference and amplifica- 
tion spectra with our wave propagation model. Inside the 
umbra of the sunspot, the magnetic field imposes wave 
propagation along field lines, and thus the hypothesis of 
vertical propagation is justified. However, in the quiet 
Sun waves can propagate in different directions and this 
approximation may not be true. Still, we were able to 
find a fit matching reasonably the phase spectrum. We 
retrieve a lower value of formation height differences and 
a higher temperature in the quiet Sun compared to the 
umbra (see Table [3]). 

The bottom panel of Fig. [15] shows the amplification 
spectrum for this case. For frequencies in the band of 2-7 
mHz, the theoretical amplification matches the observa- 
tional one. 

3.4. Comparison of the parameters of the fit with a 
model of sunspot 

From the fit to the phase and amplitude spectra of 
all the pairs of lines we obtain the temperature, the dif- 
ference between the formation heights of two lines, and 
the cooling time that best match the observational data. 
Obviously, our model is simplified and has several im- 
portant limitations. It only describes linear wave propa- 
gation in an isothermal stratified atmosphere, not taking 
into account a realistic temperature stratification of the 
sunspot atmosphere. In a real sunspot, the temperature 
varies with height, so the temperature that we obtain 
from the fit represents a mean value between two heights 
of formation. In the case of spectral lines formed at a 
similar height (i.e., the Fel lines), the temperature that 
we find should be close to the real temperature in the 
layer. However, when the formation height difference 
(Az) obtained from the fit is larger (i.e., Sil-Hel), we 
can not assign the temperature to a certain height. The 
presence of shocks in the temporal evolution of the He I 
and Call H velocities (Fig. [7j), together with the difficul- 
ties of the model to fit the amplification spectra between 
Call H and He I (Fig. fT4j) . show that the propagation 
at these heights is non-linear. Thus our determination 
of the height difference between these two layers must 
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Fig. 15. — Phase, coherence and amplification spectra between 
Fel A 3969.3 and Sil in a quiet Sun region. The format of the 
figure is the same as Fig. [9] 

be somewhat affected by the deviations from the linear 
regime. 

Despite these limitations, we plotted the deduced tem- 
perature values ov er the tempera t ure st ratification of the 
sunspot model of iMaltbv etaLl (|1986l ) (Fig. [J6]). The 
symbols in this figure mark the values retrieved from the 
fit to our observations. From the fit, we obtain the rel- 
ative difference between the formation heights of spec- 
tral lines, not the absolute values. As reference point, 
we set the height of the Si l line to z=308 km, as given 
by iBard & Carlssonl (|2QQ8[) f or the sunspot atmosphere 
model of IMaltbv et al.l (|1986| ). The formation heights of 
the other lines then follow from their relative distance to 
the Sil line (Table [3]). Table |4] lists the resulting height 
values and the corresponding temperature. 

In the case of the Fe I lines, we obtained from the low 
phase difference (top panel of Fig. [11]) that the distances 
between them are very small, but the amplification spec- 
tra (bottom panel of Fig. [TT]) or the histograms of ve- 
locities (Fig. [6]) show a certain increase of the amplitude 
with height. This amplification in spite of the small ge- 
ometrical distance indicates a low pressure scale height 
Ho and consequently a low temperature T. The height 
of the Fe I line close to the temperature minimum agrees 
with the results from the fit of the model to the phase 
and amplification spectra. 

In the case of the chromospheric signals (Call H core 
and He I line), we assign to both formation heights the 
temperature retrieved from the fit of the phase spectra 
between Call H and He I velocities. According to the ge- 
ometrical differences of Table [3j there is an uncertainty 
in the formation height of these two lines. On the one 
hand, we have set them from their Az with respect to 
the Sil line, locating the formation height of the Call H 
core at 958 km and the formation height of the He I line 
at 1208 km. On the other hand, we have obtained an 
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Fig. 16. — Temp erature stratification in the sunspot model of 
I Malt bv et al.l ((1986). Asterisks represent our best-fit values to the 
phase difference and amplification spectra of the line pairs. Dia- 
monds mark an alternative estimate of the formation height of the 
Call H core and the He I line. 

alternative height for the Ca II H core by considering its 
Az to the He I line, and subtracting it from the previ- 
ous estimate of the He I formation height. In the same 
way, we have located the formation height of the He I line 
taking into account the geometrical differences between 
the pairs Sil-Fel A 3969.3 and Fel A 3969.3-Hel. Thus, 
we retrieved a range of heights for both chromospheric 
lines, as it is shown in Table HJ The hei ght ranges we find 
are comparable to prev i ous works (e.g. , Lites et al.lll993l: 
Carlsson fc Steinl 119971 : ICenteno etatl 120091 : iBeck et all 
2009|), even if some of these articles deal with the quiet 
Sun solely. 



TABLE 4 

Formation heights and temperature for the spectral lines 



Spectral line 


z [km] 


T[k] 


Si i 


308 


4000 


Fei 3965.4 


558 


3500 


Fei 3966.0 


568 


3500 


Fei 3969.3 


588 


3500 


Can H 


958-1108 


6000 


He I 


1038-1208 


6000 



4. DISCUSSION AND CONCLUSIONS 

We have presented an analysis of the LOS velocities 
obtained from at set of spectropolarimetric data in the 
near-IR spectral region around 10830 A and the optical 
region around 3969 A in a sunspot atmosphere and its 
vicinity. From these two spectral ranges, we retrieve on 
the one hand the Doppler velocities of the photospheric 
Si I A 10827 and the chromospheric He I A 10830 line. 
On the other hand, we also sample several layers be- 
tween these two heights, using the Doppler shifts of the 
chromospheric Call H line core and the photospheric Fel 
lines from the wings of the Call H line. 

The histograms of LOS velocities show that the width 
of the velocity distribution increases with height, both in 
the sunspot and the quiet Sun atmosphere. Quiet Sun 



rms velocities are larger than those in the sunspot due 
to the higher power of oscillations in the quiet Sun. The 
growth of the amplitude of the oscillations with height 
is scaled by the pressure scale height Hq. It is smaller 
in the umbra than in the quiet Sun atmosphere, so the 
amplitude of umbral oscillations rises faster. This yields 
that the ratio of the quiet Sun and umbral rms velocities 
of the same lines decreases with the formation height of 
the lines. 

The phase difference spectra of LOS velocities between 
several pairs of lines show upward propagating waves for 
frequencies higher than 4 mHz. The power at lower fre- 
quencies does not propagate up, since waves with these 
frequencies are evanescent. The slope of the phase spec- 
tra, together with the histograms of LOS velocity, allows 
us to sort all the spectral lines used in this work by their 
formation height. 

Phase and amplification spectra were fitted to a model 
of linear wave propagation in a stratified atmosphere 
with radiative losses following the Newton cooling law. 
The model works reasonably well at layers below the 
formation height of the Call H line core, where waves 
propagate in a linear regime, while it fails in the fit of 
the amplification spectrum between the He I velocity and 
the Call H line core velocity, due to the importance of 
nonlinear it ies at these chromospheric heights. 

From the fit to the model, we retrieved the temper- 
ature, the difference in geometrical height between the 
formation heights of both spectral lines, and the radia- 
tive relaxation time. Setting the hei ght of the lower- 
most forming line (Sil) to agree with IBard fc Carlsson! 
(|2008[ h the formation height of all the lines in sunspots 
was inferred. The Fe I lines from the wings of the Ca II H 
line are formed about 250 km above the photospheric Si I 
line. The relative position of the lines is well determined 
by the observations, since the rms velocities (Fig. [6]), 
the power spectra (Fig. |8j), and the phase and amplifi- 
cation spectra (Fig. [T2]) all indicate that the Fel lines 
are formed in the upper photosphere above the Sil line. 
The temperatures obtained for the spectral lines then 
show a goo d agreement with th e temperature stratifica- 
tion of the iMaltbv et al.l ([ 19861 ) sunspot model, and the 
formation heights are coherent with previous estimates. 
It must, however, be taken into account that the esti- 
mate of the formation h eight from response functions in 
IBard fc Carlssonl ([20081 ) was performed for a static atmo- 
sphere and has an strong dependence on the atmosphere 
model employed. 

Most of the power of the photospheric lines is concen- 
trated in a prominent peak at 3.5 mHz, in the 5 minute 
band. From a comparison between the power spectra 
of the Sil and one of the Fel lines, it is interesting to 
note that the power peak is exactly at the same position, 
although the iron line forms at around 200 km above 
the silicon one. The same behavior was found for the 
chromospheric He I line and the Call H line core, which 
peak at around 6 mHz, corresponding to the 3 minute 
band. The maximum of the power spectra is not shifted 
gradually from 3.5 mHz at the photosphere to higher fre- 
quencies at larger heights, but the photospheric and chro- 
mospheric group of lines show a discontinuous behavior. 
Waves at frequencies above the cut-off increase their am- 
plitude with height faster than evanescent waves below 
the cut-off, resulting in larger power of 3-minute waves 
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at chromospheric heights. Yet, as follows from Fig. [Hb, 
the LOS velocity power at frequencies above 4 mHz mea- 
sured from of the Fe I line is higher than the one from the 
Si I line; it means that already in the upper photosphere 
the high-frequency power becomes important. This find- 
ing suggests that high-frequency waves prominent in the 
chromosphere have to be generated in the photosphere or 
below and their dominance at the chromospheric height is 
the result of their large amplitude increase. These results 
are consistent with t hose obtained from observations by 
iCenteno et al.l ([200 6) using t wo He I and Si I line s, and 
from numerical simulations by iFelipe et al.l (|2010l ) . 
From the compatibility between our observations and a 



simple wave model, we conclude that we observe a contin- 
uous field-aligned propagation of slow magneto-acoustic 
waves in the upper atmosphere of sunspot. These waves 
first reach the formation height of Si I, then the formation 
height of the Fel lines from the Call H line wing located 
in the upper photosphere, then the formation height of 
the Call H line core and finally that of the He I line. The 
propagation becomes non-linear at heights between the 
formation of Fel lines and Call H line core. 
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